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O ■ 

^ ■ ABSTRACT 

Context. Ages, metallicities, space velocities, and Galactic orbits of stars in the Solar neighbourhood are fundamental observa- 
tional constraints on models of galactic disk evolution. Understanding and minimising systematic errors and sample selection 
. biases in the data is crucial for their interpretation. 
^ ' Aims. We aim to consolidate the calibrations of uvbyP photometry into Tca, [Fe/H], distance, and age for F and G stars and 
T-H ' rediscuss the results of the Geneva-Copenhagen Survey (Nordstrom et al. 2004; GCS) in terms of the evolution of the disk. 
G^ Methods. We use recent V-K photometry, angular diameters, high-resolution spectroscopy, Hipparcos parallaxes, and extensive 
numerical simulations to re-examine and verify the temperature, metallicity, distance, and reddening calibrations for the uvhyj3 
, system. We also highlight the selection effects inherent in the apparent-magnitude limited GCS sample, 
f — ■ Results. We substantially improve the Teff and [Fe/H] calibrations for early F stars, where spectroscopic temperatures have 
■ large systematic errors. A slight offset of the GCS photometry and the non-standard helium abundance of the Hyades invalidate 
, its use for checking metallicity or age scales; however, the distances, reddenings, metallicities, and age scale for GCS field stars 
' require minor corrections only. Our recomputed ages are in excellent agreement with the independent determinations by Takeda 
^ I et al. (2007), indicating that isochrone ages can now be reliably determined. 
• ^ . Conclusions. The revised G-dwarf metallicity distribution remains incompatible with closed-box models, and the age-metallicity 
' relation for the thin disk remains almost fiat, with large and real scatter at all ages {a-i„trinsic~ 0.20 dex). Dynamical heating of 
JH , the thin disk continues throughout its life; specific in-plane dynamical effects dominate the evolution of the U and V velocities, 
. . . 1 while the W velocities remain random at all ages. When assigning thick and thin-disk membership for stars from kinematic 
criteria, parameters for the oldest stars should be used to characterise the thin disk. 

Key words. Galaxy: stellar content - Galaxy: solar neighbourhood - Galaxy: disk - Galaxy: kinematics and dynamics - Galaxy: 
evolution - Stars: fundamental parameters 



1. Introduction 

Models for the evolution of spiral galaxy disks describe 
their star formation history, nucleosynthesis, chemical en- 
richment, and dynamical evolution. Traditional parame- 
terised models yield single- valued relations for the increase 
in the (total or individual) heavy-element abundances of 
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stars at a given position in the disk, radial gradients in ele- 
mental abundances, and the kinematic heating of the local 
disk, all as functions of time. The underlying paradigm is 
the efficient mixing and recycling of interstellar gas, such 
that mean values of abundances and kinematics as func- 
tions of time and radial position in the disk describe the 
underlying physical processes adequately (see Casuso & 
Beckman I^ODil Naab & Ostriker [M^Sl or Cescutti at al. 
120071 for recent examples). Any dispersion of the observed 
values around the mean relations is assumed to be due to 
observational error. 

Under realistic conditions, however, local variations in 
the rate of chemical enrichment must have occurred (see 
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e.g. Brook et al. 120071 for such an approach). The questfon 
of interest is therefore how well these relations and their 
intrinsic scatter can be determined from the observations. 

The Milky Way is the one galaxy in which these pre- 
dictions can be tested in detail. Thus, complete, accurate 
information on the stellar content of the Solar neighbour- 
hood remains a fundamental observational constraint on 
any set of models. Much of the discussion focuses on the 
age-metallicity relation (AMR) for the Solar neighbour- 
hood, and the key questions are two-fold fTwarog I1980[ 
Carlberg et al. 119851 Meusinger, Stecklum & Reimann 
fT99T| Edvardsson et al. fT993l Rocha-Pinto et al [20001 
Feltzing, Holmberg & Hurley [200 H and Nordstrom et al. 
I2004p : (i) Does the average relation show the expected rise 
in metallicity from the formation of the (thin) disk to the 
present time? and (ii) how large is the intrinsic dispersion 
in metallicity at any given age? Given the wide-ranging 
ramifications of these questions, the observational data 
used for the test must be prepared, selected, and discussed 
with the utmost care. 

The most comprehensive recent study of nearby stars 
in the solar neighbourhood is the Geneva-Copenhagen 
Survey (Nordstrom et al. [2001 GCS in the following). 
The GCS provides metallicities, ages, kinematics, and 
Galactic orbits for a complete, magnitude-limited, all- 
sky sample of ^^14,000 F and G dwarfs brighter than 
V ~ 8.3. The basic observational data are uvbyP pho- 
tometry, Hipparcos/Tycho-2 parallaxes and proper mo- 
tions, and some 63,000 new, accurate radial velocity ob- 
servations, supplemented by earlier data. The best cali- 
brations then available were used to derive Tcff, [Fe/H], 
and distances from the photometry. The astrometry and 
radial velocities were then used to compute space motions 
and identify binaries in the sample. Finally, unbiased ages 
and error estimates were computed from a set of theo- 
retical isochrones by a sophisticated Bayesian technique 
(J0rgensen & Lindegren I2005p . and Galactic orbits were 
computed from the present positions and velocity vectors 
of the stars and a Galactic potential model. 

The GCS is large enough to yield adequate statis- 
tics for subsets of stars defined by age, metallicity, or 
abundance, and is essentially free of the kinematic and/or 
metallicity biases affecting most earlier samples. However, 
systematic selection effects may still remain in the data 
or be introduced by the calibrations used to derive astro- 
physical parameters from the observations. Accordingly, 
the purpose of the present paper is to critically re-examine 
the observational determination of Teg, [Fe/H], and age 
for F- and G-type dwarf stars in the light of the most 
recent developments in the field. We also compare our re- 
sults with those of the recent papers by Haywood (|2006( 
H06 in the following), Valenti & Fischer PUDSI VF05 in 
the following), and Takeda et al. p007|) . 

We begin the paper by briefly describing the possible 
impact of systematic and random errors in these param- 
eters on the overarching theme of spiral galaxy evolution 
in Sect. O The reader mostly interested in our new uvby(3 
calibrations in terms of effective temperature, metallic- 



ity, distance, reddening, and age computations will find 
these discussed in Sect. [3] - |7l The reader primarily in- 
terested in the evolution of the Galactic disk may skip 
directly to Sect. [51 where we discuss the end-to-end sim- 
ulations we have used to verify the robustness of the re- 
sults reached with the new calibrations. We rediscuss the 
"G dwarf problem" , age-metallicity, and age- velocity re- 
lations for the Solar neighbourhood in Sect. [9] - [Tl] and 
compare the thick and thin disks in Sect. [T^l Finally, we 
summarise our findings and conclusions in Sect. [T3l 

2. Astrophysical implications of calibration errors 

The key parameters to be compared with models are the 
masses (i.e. main-sequence lifetimes), ages, heavy-element 
abundances, spatial distributions, and space velocities or 
Galactic orbits for sufficiently large samples of stars that 
are representative of the general population of disk stars. 
However, except for positions and velocities, these param- 
eters cannot be determined by direct observation. 

Typical observational data are multiband colour in- 
dices and perhaps spectra, from which Tcs, [Fe/H], and 
absolute magnitude or log g are derived, using theoreti- 
cal or empirical calibrations. From Tes, [Fe/H], and My, 
in turn, the age and mass of each star can be derived by 
comparison with stellar models. 

This is straightforward in theory. In practice, it is 
highly non-trivial because the problem is very non-linear, 
and significant uncertainties in the calibrations remain. 
In particular, uncertainties in the theory of stellar atmo- 
spheres continue to play a major role, both for the pre- 
dicted effective temperatures of stellar models and for the 
transformations between Tog, [Fe/H], and log g and the 
observed colour indices. Because of the strong correlations 
between these parameters, calibration errors in one may 
bias the determination of another and lead to spurious 
correlations between the derived quantities, and incorrect 
conclusions about the evolution of the Galactic disk. 

As one example, assume that the adopted calibration 
yields too high T^s for the hotter stars (as indeed we find 
below). [Fe/H] as determined from equivalent widths of 
spectral lines will then be overestimated. The hotter T^g 
corresponds to a lower age, the higher [Fe/H] to cooler 
models and thus an even lower derived age for the observed 
star. A spurious age-metallicity relation results. 

As another example, assume that the T^g scale of the 
models is correct for the solar abundance, but too hot 
at lower [Fe/H], as was found in GCS. This will lead to 
overestimated ages for the metal-poor stars, i.e. again a 
spurious age-metallicity relation. A similar - or additional 
- error is introduced if the enhanced [a/Fe] ratio of metal- 
poor stars (e.g. Edvardsson et al. I1993P is ignored; if the 
heavy-element abundance parameter Z of the models is 
assumed to scale simply as [Fe/H], too hot models are 
selected, and the resulting age is overestimated again. 

Finally, the importance of a clear understanding of the 
definition of the stellar sample cannot be over-emphasised: 
Taking the AMR as the prototype again, the choice of stel- 
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lar sample can determine the outcome even before a sin- 
gle observation is made. E.g., the F-dwarf sample studied 
by Edvardsson et al. (|1993p excluded any old, metal-rich 
stars that might exist and a priori decided the shape of 
their AMR - a fact emphasised in the paper, but largely 
ignored in later references. More complete and accurate 
observational data are always useful, but temperature cal- 
ibrations, the choice and verification of stellar models, and 
the methods used to compute ages and their uncertainties 
are far more urgent issues at present. 

3. Temperature calibration 

Teft is the most critical parameter in the determination of 
isochrone ages, but also affects the spectroscopic metal- 
licity determinations. Two qualitatively different methods 
to determine TeS exist, based on measurements of spec- 
tral lines or on colour indices. When using spectra, the 
determination of Toff is usually based on the excitation 
balance of iron, but depth ratios of sets of spectral lines 
have also been used. When using photometry, the deter- 
mination of Toff is usually tied to the infrared flux of the 
stars as estimated from models (the IRFM technique). 

For a fundamental test, we return to the basic defini- 
tion of Toff: /boi = ^cr'Toff' where /boi is the bolometric 
flux and the angular diameter of the star. The prac- 
tical problem in applying the formula is that solar-type 
main-sequence stars have very small angular diameters 
that are difficult to measure accurately. The situation is 
quickly improving, however, as new interferometers enter 
operation, notably the ESO VLTI. Kervella et al. (|2004|) 
summarise the situation and give diameters for 20 A-M 
main-sequence stars and 8 A-KO sub-giants. The angular 
diameters are of excellent quality, especially those from 
the VLTI, which have errors down to 1%. 

Ramirez & Melendez (|2005aP combined these diam- 
eters with bolometric flux measurements to derive Toff 
directly for 10 dwarfs and 2 sub-giants. Propagating the 
errors in the diameters and fluxes through to the tempera- 
tures, the mean error is 57K. Ramirez & Melendez (j2005ap 
also give (IRFM) estimates for these stars, with a mean 
difllerence < Tlf^ - T^^' >= lOK, dispersion 98K, and 
mean error 28K. 

As regards spectroscopic Teff determinations, Santos 
et al. (|2005l and references therein) give results for stars 
both with and without detected planetary companions, 
including 7 stars with measured angular diameters. The 
mean difference is < T^^^*""" - T^^^' >= 92K, with a dis- 
persion of 91K and mean error 34K. This is in agreement 
with Santos et al. (i2005p themselves, who flnd their Toff 
scale to be 139K hotter than that by Alonso et al. p996p . 
based on the IRFM and calibrated to B-V and [Fe/H]. 

Another valuable comparison is with the recent paper 
by VF05. They give spectroscopic data for ^1000 stars 
on the Keck/Lick/ AAT planet search program, based on 
fits of synthetic spectra to their high-resolution data and 
correcting their zero-points from observations of Vesta as 



a proxy for the Sun. For the 8 stars with measured diam- 
eters, the mean difference is < Toff^'" -T„^# 64K, with 
a dispersion of 124K and mean error 44K. 

The reasons for the failure of spectroscopic Tcs deter- 
minations based on the excitation equilibrium in ID static 
LTE models are given in e.g. Asplund (1 2005]) . Basically, 
because real stars are spherical, hydrodynamical systems 
with lines formed in NLTE, the three fundamental as- 
sumptions underlying traditional model atmospheres are 
inadequate and lead to biased results. The correction pro- 
cedure used by VF05 seems to have eliminated at least 
part of this bias. 

In the following, we review the GCS and other recent 
temperature scales and derive an improved calibration for 
the uvby(3 system. 

3.1. Teff in the GCS 

The marked offset of the Santos et al. (|2005p Teff values 
is confirmed for the 160 stars common to Santos et al. 
and the GCS, which employed temperatures based on the 
Alonso et al. p996p calibration of b-y, mi and ci to the 
IRFM scale. The mean difference is < T^f"*"'' - T^GC >^ 
127K (dispersion 84K, mean error 7K). For the spectral 
type range in question, this shows that the GCS Toff scale 
is in excellent agreement in the mean with temperatures 
derived directly from the definition of Toff. 

However, the Alonso et al. p996p calibration of b-y, 
TOi and ci may be valid only in the ranges in Teff within 
which enough calibration stars exist. Fig. 11a of Alonso 
et al. (|1996p shows the fitted temperatures as a function 
of b-y and makes clear that there is a marked dearth of 
calibration stars blueward of b-y « 0.3 or Toff « 6500K; 
a similar lack is seen for the reddest stars. Alonso et al. 
(|1996p give a dispersion for their relation of cr(0off) = 0.019 
(llOK), while Ramirez & Melendez fj2005b;) give 87K for 
their new b-y relation, albeit at the price of not being valid 
over the whole parameter range of the GCS stars. 

The potential problem in the Alonso et al. (|1996p uvby 
calibration for blue and red stars is clearly demonstrated 
by a comparison with other calibrations of Stromgren pho- 
tometry. They all agree rather well between 0.3<b-y<0.6, 
but diverge outside these hmits, see e.g. Fig. 24 in Clem 
et al. ([2004]) . 

The situation is unfortunately similar for stars with 
direct Toff determinations, where Kervella et al. (|2004p 
have only one star between 6000 and 8500K (Procyon at 
6500K), compared to 5 stars between 8500 and lOOOOK 
and 11 between 5000 and 6000K. 

3.2. Calibrations based on V-K 

Alternatively, temperatures based on V-K can be used, 
because 2MASS photometry is available and well suited to 
our colour and magnitude ranges: Of the 16682 GCS stars, 
16139 have 2MASS K-magnitudes of the highest quality 
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0,2 0.3 0.4 0,5 0.6 

b-y 

Fig. 1. Differences between Toff from the Alonso et al. 
P996p b-y calibration, used in the GCS, and Teg as de- 
rived from the V-K calibration by di Benedetto (|1998p . 

class "A" . We combine these with the V magnitude from 
the GCS uvhy[3 photometry. 

In contrast to calibrations based on only visible 
colours, Tcff estimates based on V-K are remarkably ro- 
bust; the small differences are probably mostly related to 
the different existing definitions of the K-band. As an ex- 
ample, the V-K calibrations of Alonso p996p . Ramirez & 
Melendez H20Q5bp . and di Benedetto (|T998l) agree pairwise 
to within ~20K for the stars in the GCS. In contrast, the 
Alonso V-K and h-y calibrations differ by 180K rms! 

Today, the K-band system of choice is that defined by 
the 2MASS catalogue, which also gives K magnitudes of 
high quality for most of the stars in the GCS. We have used 
the calibration of di Benedetto (|1998p . which is in excellent 
agreement with the IRFM temperature scale of Ramirez & 
Melendez (|2005ap and is valid for the whole colour range of 
the GCS. In order to transform the Johnson V-K colours 
used by di Benedetto to the Ks system used by 2MASS, 
we use the relation: 

{V - K)j = 1.007[(y -Ks- 0.044) - 0.01] 

Fig.[T]shows the difference between the T^s values from 
the original GCS catalogue and those derived from the 
2MASS photometry and the corresponding calibration by 
di Benedetto (|1998p . The systematic differences are quite 
marked, especially for the hotter stars. 

The trend seen in Fig. [1] is qualitatively similar to that 
found by H06 (Fig. 7), also by comparison with values 
from di Benedetto (|1998p . but is better determined here, 
due to our much larger sample and because V-K is a much 
more reliable temperature indicator than B-V. 

3.3. A new b-y - T^s calibration 

Because of the larger number of suitable calibration stars 
available, the V-K calibration has less systematic error. 
However, the observational accuracy of the b-y index is 
substantially better than that of V-K. We have therefore 
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Fig. 2. Differences between the Toff of GCS stars as de- 
rived from our new b-y calibration and from the V-K cal- 
ibration by di Benedetto (|1998p . 



derived a new Toff calibration for b-y, based on the V-K 
temperature scale set by di Benedetto (1998). The best 
result is achieved when the sample is divided into three 
temperature ranges. In the blue range (0.20<&-?/<0.33) 
the relation is: 

6(off = 0.649 + 0.095(5 - y) + 0.034[Fe/H] -f 1.207(6 - y)^ 

- 0.005[Fc/H]^ - 0.181(6 - y)[Fc/H] 

where 6'off — 5040K/TefT. The dispersion of the fit is 
<^{&cs) — 0.007 when 2.5-a outhers are removed, corre- 
sponding to 60K or 0.004 in logTgff for the mean log Toff. 
In the middle range (0.33<6-?/<0.50), the relation is: 

Oes = 0.754 - 0.365(6 - y) - 0.001[Fe/H] + 1.635(6 - yf 

- 0.011[Fc/H]^ - 0.091(6 -y)[Fc/H] 

The dispersion of the fit is CT(6'eft) — 0.009 when 2.5-ct 
outliers are removed, corresponding to 57K or again 0.004 
in log Toff for the mean log Toff. 
In the red range (0.50<6-?/<0.60), the relation is: 
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(9eff = 0.255 + 1.656(6 - y) + 0.018[Fe/H] - 0.397(6 - y)^ 
-0.011[Fc/H]^ -0.101(6 -?/)[Fc/H] 

The dispersion of the fit is a{9cs) = 0.012 when 2.5- 
a outliers are removed, corresponding to 55K or 0.005 in 
logTeff for the mean log Toff 

We note that this calibration yields 5777K for the 
Sun when using (b-y)Q= 0.403 from Holmberg, Flynn & 
Portinari (|2006p . although this was not forced on the fit. 
The two main calibrations are also very well connected: At 
their common colour, (b-y)— 0.33, the difference in tem- 
perature is 2, 6, and 7K at [Fe/H]= 0, -0.50, and -1.00. 
Due to the small number of very red calibration stars, the 
difference in Tes grows to -19, 18, and 52K at [Fe/H]= 0, 
-0.50, and -1.00 at (h-y)= 0.50 - stiU very smaU. 

Fig.[2]demonstrates the agreement of the resulting Tcff 
values with those derived from V-K and the di Benedetto 
(figgg)) calibration. 

Finally, we compare our new photometric tempera- 
tures to the spectroscopic results by VF05. Fig. [3] shows 
the temperature differences for the 697 single stars in com- 
mon with the GCS. The mean difference is similar to that 
from the stars with angular diameters: < Trff^^— T^J^g ^ >— 
43K, with a dispersion of 91K. 



0.4 
b-y 



Fig. 3. Differences between Tcff from our new b-y calibra- 
tion and Teff as derived from spectrum fitting by VF05. 



4. Metallicity calibration 

[Fe/H] is important, both in itself as a diagnostic of 
the chemical evolution of the disk, and because it en- 
ters into the determination of stellar ages from theoretical 
isochrones. However, a comparison of the GCS photomet- 
ric metallicities with the several high-quality spectroscopic 
studies available today (Fig. |4]) shows significant differ- 
ences, even when considering only spectroscopic studies 
using a photometric temperature scale. 

For the four major spectroscopic studies (in the order 
Edvardsson et al. [T9931 Chen et al. [20001 Reddy et al. 



■■■mm; 
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[Fe/H] 



Fig. 4. Modern spectroscopic [Fe/H] determinations vs. 
the photometric [Me/H] of the GCS catalogue. 



120031 and AUendc Prieto et al. l2004p . the mean differences 
< [Fe/H] — [Me/H] > are, for the F-calibration range: 
0.064, 0.038, 0.114 and 0.013; and for the G-calibration 
range: 0.019, -0.038, 0.022 and -0.044. Only AUende Prieto 
et al. (|2004p have stars in the range of the "K" calibration 
(as defined in GCS), with < [Fe/H] - [Me/H] >= -0.005. 

As regards Santos et al. (|2005p . who use a spectro- 
scopic Teff scale, we find < [Fe/H] - [Me/H] >= 0.162 
for stars in F range, < [Fe/H] - [Me/H] >= 0.082 for 
stars in the G range, and < [Fe/H] - [Me/H] >= 0.022 
for their stars in the K range. We conclude that metal- 
licities derived using the spectroscopic Toff scale must be 
reduced by about 0.10, 0.09 and 0.03 dex for the F, G and 
K temperature ranges. This is in fair agreement with the 
estimated effect on [Fe/H] of a temperature shift of about 
lOOK (see e.g. Fehzing & Gustafsson [19981) ■ 

We conclude that, when comparing metallicity scales 
from different sources, it is crucial to verify how these 
scales are established. It is especially important to check 
which Toff scale is used, because an erroneous Toff scale 
can introduce large biases in the abundances, especially 
for solar or somewhat hotter temperatures. 

4.1. The original GCS calibration 

The determination of accurate metallicities for F and G 
stars is one of the strengths of the Stromgren uvbyf3 sys- 
tem. Among the then available uvby(3 calibrations, the 
GCS used that by Schuster & Nissen (|1989p for the ma- 
jority of the stars. 

However, this calibration was found to give substantial 
systematic errors for the very reddest G and K dwarfs 
{b — y > 0.46), where very few spectroscopic calibrators 
were available at that time. In the GCS we therefore 
derived a new relation, based on a sample of 72 dwarf 
stars in the colour range 0.44 < b — y < 0.59, using the 
same terms in the calibration equation as the Schuster & 
Nissen (|1989p G-star calibration. The resulting equation 
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[Fe/H] = -2.06 
+ 73.50m2(6 - 



f 24.56mi - 31.61mf - 53.64toi(6 - y) 
y) + [26.34TO1 - 0.46ci - 17.76m2]ci 



The metallicities from this cahbration are compared 
to the new spectroscopic reference values in Fig. [S] The 
dispersion around the (zero) mean is 0.12 dex. 

For the --600 stars in the interval 0.44<&-?/< 0.46, the 
new calibration agreed with that by the Schuster & Nissen 
(fT989)) to within 0.00 dex in the mean (s.d. 0.12). 

About 2400 GCS stars with high T^s and low log g 
were outside the range covered by the Schuster & Nissen 
(|1989p calibration. For these stars, the calibration of j3 and 
mi by Edvardsson et al. (|1993p was used when valid. For 
the stars in common, the two calibrations agree very well 
(mean difference of 0.00 dex, dispersion only 0.05). 

For stars outside the limits of both calibrations, we 
derived a new relation, using the same terms in the 
equation as the Schuster & Nissen (|1989p calibration for 
F stars. From 342 stars in the ranges: 0.18 <b — y< 0.38, 
0.07 < mi < 0.26, 0.21 < a < 0.86, and - 1.5 < [Fe/H] < 
0.8, we found the following calibration equation: 



[Fe/H] = 9.60-61.16mi 

,2 



81.25mi{b~y) 
224.65m^(5 - y) - 153.18mi(6 - yf 
[12.23 - 90.23mi + 38.70(5 - y)] log(mi 



C3), 



where 
of the 



C3 = 0.45- 3.98(6 -j/) + 5.08(6 -j/)2. The dispersion 
fit is 0.10 dex. 




[Fe/H] 



Fig. 5. GCS photometric metallicities ([Me/H]) vs. the 
spectroscopic [Fe/H] values used to establish the calibra- 
tions. Open circles depict the cool stars, dots the hot stars 
(see text). 



Fig.Ocompares the photometric [Fe/H] from this cali- 
bration with the spectroscopic values; for the stars in com- 
mon, they again agree very well those from the Schuster 
& Nissen (|1989|) calibration (mean difference 0.02 dex, a 
only 0.04). 




b-y 



b-y 



Fig. 6. Left: Differences between the GCS metallicities for 
all stars and those now derived from the GCS blue calibra- 
tion (the sharp line at A [Fe/H] = corresponds to those 
GCS stars that already used this calibration). Right: The 
difference between metallicities from the new calibration 
derived in this paper and the original GCS metallicities, 
based on the Schuster & Nissen ()1989|) cahbration. 



4.2. An improved metallicity calibration 

In order to define an improved photometric metallicity 
calibration that is valid for all GCS stars, we select only 
spectroscopic investigations using a photometric temper- 
ature scale (unlike the compilation of literature values by 
Cayrel de Strobel et al. ([MOT]) used by H06). This yields 
a total sample of 573 stars from Edvardsson et al. (|1993p . 
Chen et al. [20001 Reddy et al. [20031 AUende Prieto et al. 
(PUni|) . and Feltzing & Gustafsson (fTM5| . They aU have a 
very consistent metallicity scale: Comparing stars in com- 
mon yields differences of of 0.02 dex or smaller in the 
mean, with a dispersion of ^0.07 dex. The 573 stars span 
the range 0.24 < b-y < 0.63, 0.10 < mi < 0.70, 0.17 
< ci < 0.53, and -1.0< [Fe/H] < 0.37. 

A new fit of the uvby indices to the spectroscopic 
[Fe/H] values from this sample was performed, using a 
calibration equation containing all possible combinations 
of 6 — mi, and ci to third order. 

The resulting calibration equation is: 

[Fe/H] = -2.19 - 1.02(6 - y) + 7.34mi - 0.27ci 
+ 5.86(6 - y)2 - 43.74mf - O.Mcf + 25.03(6 - y)mi 
+ 5.29(6 - y)ci + 25.95miCi - 31.10(6 - y)^ + 46.19m? 



-3.86cf-H4.54(6-y)^m: 
'18.36m?ci-9.99c?(6- 



+ 19.31(6-y)2ci - 17.46m2(6-y) 
- y) + 6.60c^mi — 59.65(6 — y)mici 



This new calibration is used for all stars with 0.30 < b-y < 
0.46. Together with the blue and red calibrations already 
derived in the GCS, this completes the final metallicity 
calibration of the present paper. All results discussed in 
the following are based on it. 

Fig. |S| shows the differences between the GCS [Fe/H] 
values and those derived with the new calibration; note 
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Fig. 7. Same as Fig. |4l but using the new metallicity cal- 
ibration derived in this paper. 
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Fig. 8. Spectroscopic [Fe/H] data from VF05 vs. our new 
photometric [Me/H] for the GCS stars. 



that, in the left panel of Fig. [HI the heavy sequence at 
zero difference is formed by the stars that already then 
used the blue calibration. 

The consistency of the photometric metallicities from 
this calibration with the spectroscopic reference values 
is shown in Fig. [T] The dispersion around the mean is 
0.07 dex - the same as between two different spectro- 
scopic measurements. The calibration derived above is 
also in good accordance with the two relations derived 
in the GCS. Compared to the blue relation, the mean 
difference is 0.02 (dispersion 0.04) in the common range 
0.3<b-y<0.32. Compared to the red relation, the mean 
difference is 0.00 (dispersion 0.09) in the common range 
0.44<6-y<0.50. 

Finally, Fig. [8] compares our new GCS metallicities 
with the spectroscopic [Fe/H] values from VF05. The 
mean difference VF05-GCS is 0.08 dex, with a standard 
deviation of O.IO dex. The higher metallicities derived by 
VF05 are understandable, given the generally higher tem- 
peratures they adopt. 



Fig. 9. Top left: The GCS metallicity calibration applied 
to the Crawford & Perry (|1966p uvhyp photometry of the 
Hyades. Top right: Using the new calibration from this pa- 
per instead. Bottom left: The GCS metallicity calibration 
applied to the Crawford & Barnes (1969) uvhy(3 photom- 
etry of the Coma cluster. Bottom right: The Coma stars, 
using the new calibration. 



4.3. The Hyades and Coma clusters 

One of the two open clusters included in the GCS, the 
Hyades, has been used by e.g. HOG to assess the zero- 
point of the GCS metallicity calibration. A comparison 
with the original uvby(3 photometry of the Hyades and 
Coma clusters by Crawford & Perry l|1966p and Crawford 
& Barnes (1969) shows, however, that the photometry of 
the cluster stars that was listed in the GCS is not on the 
same system as the rest of the catalogue, perhaps because 
they were observed at high air mass from Chile. While 
unfortunate, this offset means that the Hyades data will 
give misleading results if applied to the metallicities (and 
thus ages) of the whole GCS catalogue. 

The result of using the photometry from the origi- 
nal sources combined with the earlier GCS calibrations 
is shown in Fig. [5] For the Hyades, we then find [Fe/H] — 
0.10 ± 0.10; for Coma, [Fe/H] = -0.02 ± 0.08 (mean 
and dispersion). The derived [Fe/H] for the Hyades cor- 
responds well to the standard value, as expected because 
it is mostly based on the Schuster & Nisssen calibration, 
which used the Hyades as a main anchor point. 

Using the new calibration instead yields [Fc/H] = 
0.06 ± 0.07 for the Hyades and [Fe/H] = -0.09 ± 0.06 
for Coma. The dispersion among the cluster stars is no 
larger than in the field star calibration, but there is an 
offset of about 0.07 dex. 
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A different photometric [Fe/H] for the Hyades is to be 
expected, considering the non-standard He/Fe ratio of this 
cluster (Vandenberg & Clem [20051) . but Coma has a stan- 
dard chemical mix, and the mean photometric metallicity 
is very near the accepted spectroscopic value. 

In summary, we conclude that it is inadvisable to use 
the Hyades for general calibrations of photometric metal- 
licities or comparing metallicity scales. The same is true 
for the ages, especially for the unevolved solar-type stars 
for which isochrone ages are bound to be very uncertain 
by any method. 

Finally, we note that the revised data in Table 1 for 
Hyades and Coma stars are based on the standard pho- 
tometry and can be used with confidence. 

4.4. Solar analogs as a test of the calibrations 

Solar analogs offer yet another test of the temperature and 
metallicity determinations. Such stars can be identified 
through either photometric or spectroscopic resemblance 
to the Sun - preferably both. 18 Sco (HD 146233) is one of 
the best known examples, with b-y— 0.404, My— 4.77, as 
compared to the value b-y— 0.403 for the Sun by Holmberg 
et al. (j2006p . determined by comparison with a carefully 
selected sample of stars similar to the Sun. 

[Fe/H] for 18 Sco is close to zero with both the old 
(-1-0.03) and new (-0.02) metallicity calibrations, whereas 
Teff= 5766K as derived from b-y and the new calibration 
is much closer to that of the Sun than the CCS value of 
5689K, in good agreement with the general trend of Fig.[T] 

5. Absolute magnitude/distance calibration 

In order to determine an improved absolute magni- 
tude/distance calibration for uvby(3 photometry, we se- 
lected a sample of 2451 stars with absolute magnitudes 
better than 0.10 mag from Hipparcos and no indication 
of binarity in the CCS. A relation using all combinations 
of b-y, mi and Ci up to third order was then fitted to the 
absolute magnitudes of these stars. 

The difference between the calibration data and the 
fitted relation is shown in Fig. [TOl The dispersion of the 
fit is 0.24 mag, but there is a clear excess of stars with 
large differences, probably due to still-undetected binaries. 
Repeating the fit after removing stars more than 2.5 a 
from the mean relation reduces the dispersion to 0.16 mag. 
The stars span the range: 0.20 < b-y < 0.60, 0.09 < mi < 
0.64, 0.16 < ci < 0.80, and -0.84 < My < 7.24. 

The resulting calibration equation, valid within the 
above parameter ranges, is: 

My = 2.99 + 2.68(6-y) + 10.08mi + 3.14ci + 15.27(6-y)2 

- 17.73m? -I- 26.28c? - 12.73(5 - y)mi - 20.79(5 - y)ci 

- 38.86mici - 17.47(6 - yf - 2A.Uml - 8.16c? 

- 56.98(5 - yfmi + 27.86(6 - yfci + 103.32m?(6 - y) 

- 161.55m?ci - 131.37c?(6 - y) + 29.51c?mi 
+ 257.37(5- ?/)miCi 



cr(My) corresponds to distance errors of only 11% and 
7% (with and without 2.5-(j outliers). Note that this is the 
total dispersion, including also the error in the Hipparcos 
distances. A further check of the new calibration based 
on the wide physical binaries in the sample confirms this 
estimate (see Sect. I7.2|) . 




0.2 0.3 0.4 0.5 0.6 0.2 0.3 0.4 0.5 0.6 



b-y b-y 



Fig. 10. Differences in My as determined from uvby(3 pho- 
tometry and from Hipparcos distances. Left: Original GCS 
photometric My. Right: My from our improved calibra- 
tion. 

This relation is a clear improvement over the photo- 
metric distance calibration used in the GCS, which had 
an uncertainty of 13% (0.28 mag). It can also be used 
to identify further binaries in the GCS catalogue, which 
only fiagged stars with a difference between the Hipparcos 
and photometric My larger than 3 tr, i.e. larger than 0.84 
mag from the photometry alone. We have not, however, 
recomputed the statistics on binaries in the GCS, as other 
sources of uncertainty remain important. 

As already expected from Fig. 10 in the GCS, there 
is no systematic difference between the GCS distances 
and those computed with the above formula (see Fig. 
[TOl) . However, the uncertainty of the individual distances 
- which dominates that of the space motions - is reduced 
enough that we have recomputed the distances and UVW 
velocity components for all the GCS stars, correcting a 
minor error in the GCS space velocities at the same time. 

6. Interstellar reddening 

In the GCS, interstellar reddenings were derived from the 
intrinsic colour calibration by Olsen (|1988p . which yields 
reddening estimates with a stated precision of 0.009 mag. 
In the GCS there are 827 stars within 40 pc which have 
E(b-y) estimates. Fig. [TT] (left) shows the histogram of 
these E(b-y) values; a best-fit pure Gaussian shows only a 
slight zero-point shift of 0.0025 and a dispersion of 0.0105 
mag. 

The distribution is quite symmetric, showing that it is 
dominated by observational errors and that real reddening 
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E(b-y) 



Fig. 11. Distribution of E(b-y) values from the GCS. Left: 
Stars within 40 pc (solid) and a Gaussian fit (dashed) for 
a mean reddening of 0.0025 mag and dispersion 0.0105 
mag. Right: Stars with distances 40 - 70 pc and the same 
Gaussian, showing a significant excess of slightly reddened 
stars. 

within 40 pc is negligible. This is no longer true for the 
distance range from 40 to 70 pc, where the distribution 
is markedly skewed, with a clear excess of true positive 
extinctions (see Fig. [Til right). Thus, contrary to the as- 
sertion by H06, this distance range contains a clear excess 
of stars with significant positive reddenings; the mean red- 
dening is 0.0048 mag for this sample. 

We conclude that the intrinsic colour calibration by 
Olsen p988p remains valid with sufficient accuracy for the 
GCS sample, so we do not recompute E(h-y) and (b-yjo. 

7. Stellar ages 

7.1. Review of stellar models 

The determination of ages for individual field stars in 
the range 1-10 Gyr is a field with rich opportunities for 
large systematic as well as random errors. Some of these 
arise from the calibrations of the observational data, as 
reviewed in the preceding chapters. Others are due to the 
various approximations used in stellar model calculations, 
especially as regards the treatment of convective core over- 
shooting in this regime of small to vanishing convective 
cores, and of the stellar atmospheres that are used to com- 
pute the observable properties of a model, notably log g 
and Teff. 

The density and sampling of the grid of models and/or 
isochrones are also important in this range, where a major 
change in isochrone morphology near the turnoff occurs for 
ages of several Gyr. Other errors again may result from the 
techniques used to interpolate in the published isochrones 
and compute masses, ages, and error estimates from the 
observed data. 

In the GCS, great effort was devoted to the deriva- 
tion of ages and masses and their uncertainties for the 
stars of the survey, using the sophisticated interpola- 
tion method and Bayesian computational techniques of 
J0rgensen & Lindegren (j2005p and accounting for the av- 
erage a-element enhancement of metal-poor stars. Exactly 
the same method was adopted in the very recent paper 
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3.74 



3.72 



3.7 



3.68 



log T., 



Fig. 12. Solar-metallicity 8-Gyr isochrones in the 
logTcff — My diagram. Thick lines: Padova (full); Victoria- 
Regina (dashed); Yale-Yonsei (dotted). Thin lines: Geneva 
models: basic (dashed) and "best" low-mass (dotted). 



by Takeda et al. (|2007L The theoretical isochrones used 
in the GCS computations were the latest set from the 
Padova group (Girardi et al. IMM Salasnich et al. lMK?]) . 
while Takeda et al. (|2007p computed an extensive set of 
new tracks using the Yale-Yonsei stellar evolution code 
(Demarque et al. l2004p . 

In order to assess the effect of different model prescrip- 
tions, we have compared the Padova isochrones to both the 
Yale-Yonsei (Demarque et al. i^0Q4j) and Victoria-Regina 
models (VandenBerg et al. l2006j) as well as to two differ- 
ent Geneva model series (the "basic" and "best low-mass" 
models, Lejeune & Schaerer l200ip . 

Fig. [13 shows that the differences between the vari- 
ous models are in fact quite small. The difference in age 
for a star located on the three isochrones in the middle, 
if derived instead from the two outlying ones, is below 1 
Gyr at the turnoff and about 2 Gyr at Mv= 5, respec- 
tively. For the 5 Gyr isochrones, the difference is about 1 
Gyr at both points; the larger percentage age spread at 
the turnoff is related to the slightly different treatment 
of convective overshooting in the different models. The 
situation is similar at e.g. [Fe/H]= -0.71, the differences 
between the isochrones being even somewhat smaller. 

Accordingly, in the following we use the same Padova 
isochrones and a-enhancement corrections as in the GCS 
itself and do not repeat the full age computations with 
all the different sets of models. This will highlight the ef- 
fect of changing the calibrations as discussed above, while 
comparing with the ages computed by VF05 and Takeda 
et al. (|2007p will illustrate the effects of photometric vs. 
spectroscopic input parameter determinations. 
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Fig. 13. Comparison between tlie observed single GCS stars with [Fe/H]= +0.23±0.05, -0.40±0.05, -0.71±0.20, and 
-1.31 ±0.30, after allowing for a-enhancement, and ZAMS Padova isochrones with (solid) and without (dotted) tem- 
perature corrections. 



J. Holmberg et al.: The Geneva-Copenhagen Survey of the Solar neighbourhood II 



11 




Fig. 14. Left: The observed single GCS stars with [Fe/H]= 0.00 ± 0.03 and the model ZAMS with the new temperature 
corrections (soHd) and without (dotted). Right: Same, for the simulated catalogue (see Sect. 18. ll and Fig. [20] on the 
sample differences). 



7.1.1. The model temperature scale 

In the GCS, close attention was paid to verifying the 
consistency between the computed lower main sequences 
and the observed unevolved stars. An appreciable offset 
was found, and substantial, metallicity-dependent, nega- 
tive temperature corrections were applied to the Padova 
models to avoid deriving spuriously large ages for the low- 
mass stars. No similar offset seems to have been found 
necessary by Takeda et al. (|2007p . perhaps because the 
Yale-Yonsei models (Demarque et al. I2004p themselves ap- 
pear to be somewhat cooler than the Padova models (see 
Fig. [T2)) . The offset is also masked by the hotter temper- 
atures and higher metallicities they adopt for their stars. 

With the new temperature and metallicity calibrations 
derived in this paper, the temperature corrections applied 
to the Padova isochroncs in the GCS also need revision. 
We find corrections of 5 log Toff = -0.005 for [Fe/H] > 0, 
-0.01 at [Fe/H]= -0.40, -0.02 at [Fe/H]= -0.70, and -0.025 
at [Fe/H] < -1.30 to be needed for the models to match 
the observed ZAMS; see Figs. fT4l and fT3l 

As a further check. Fig. [T^ compares the temperature- 
corrected solar-metallicity model main sequences with the 



distribution of stars in the observed catalogue as well as in 
a simulated sample (see Sect. 18.1] ) with errors in logTgff, 
[Fe/H], and My as observed, both in the range [Fe/H] = 
0.00±0.03. Note that the main purpose of the figure is to 
illustrate the scatter of stars around the ZAMS due to 
measurement errors; the model has not been designed to 
account for the effect of the age-dependent sample volume 
of the GCS on the distribution of stars in the HR diagram 
(cf. Sect. EH). 

Another way to check the temperature corrections is to 
consider active stars with low chromospheric ages. Fig. [T5] 
compares stars with -0.1<[Fe/H]<0.1 and chromospheric 
ages below 1 Gyr from the catalogue of Wright et al. (|2004p 
to a solar-metallicity 1-Gyr isochrone, corrected as above. 
As seen, the agreement is again quite satisfactory. 

Finally, a test can be made using single Hyades stars 
with photometry from Crawford & Perry (1966). Fig. [T6l 
compares these stars to a 0.7-Gyr corrected isochrone and 
shows the estimated ages and their uncertainties. The 
scale offset in the ages is due to the use of a standard 
Y/Z ratio in the models, which is incorrect for the Hyades 
(Vandenberg & Clem [20031). 
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Fig. 16. Left: Single Hyades stars compared to a corrected 0.7-Gyr isochrone for [Fe/H]= +0.15. Middle: Estimated 
ages and uncertainties; the scale offset is due to the non-standard Y/Z ratio of the Hyades. Right: Age-metallicity 
relation for the Hyades computed as for the whole catalogue, including only stars with age errors below 25%. 
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Fig. 15. GCS stars with chromospheric ages below 1 Gyr 
and —0.1 < [Fe/H] < 0.1, compared to a corrected 1-Gyr 
solar-metallicity isochrone. 



7.2. Consistency checks with wide binaries 

Wide physical binaries as confirmed by proper motions 
and radial velocities provide another check of the accu- 
racy of our metallicity, distance, and age determinations. 
We have selected 18 such pairs from the GCS, with sepa- 
rate measurements of all parameters and no indication of 
further multiplicity. The pairs have a wide distribution in 
[Fe/H], from -0.4 to -1-0.1, and ages in the range 0-10 Gyr. 
Fig. [T71 shows a typical example. 

With the GCS calibrations, the rms difference in 
[Fe/H] between the two components is O.lldex; using the 
new calibrations reduces it to 0.06 dex, in excellent agree- 
ment with our other error estimates. Similarly, using the 
new photometric distance calibration (Sect. O instead of 
the older ones used in the GCS reduces the average dif- 
ference in distance between the binary components from 
16±3% to 12±3% (s.e. of mean). 



The accuracy of the age determination can be esti- 
mated from the ratio of the age difference between the 
two binary components, AAge, and the combined la un- 
certainty of the two ages, a Age. For the GCS ages we 
find < AAge/aAge >— 0.87. Recomputing the ages and 
age uncertainties for the binary components, using the 
new calibrations and (lower) estimated errors for [Fe/H], 
Toff from b-y, and photometric My as well as the newly 
corrected isochrones, we find < AAge/aAge >= 0.76. In 
other words, the age estimates become more consistent as 
well as more precise with the new calibrations. 




Age (Gyr) 



Fig. 17. My, [Fe/H], and ages for the components of the 
physical pair HD 130140/1. The isochrone shown is for the 
best mean age of 9 Gyr. AAge/aAge— 0.74 for this pair. 



7.3. New ages vs. the GCS ages 

Fig. [18] compares the GCS ages, computed with the old 
[Fe/H] and Tea calibrations, with the ages computed using 
the new calibrations and model corrections from this pa- 
per. Overall, the differences are insignificant, much smaller 
than the estimated individual uncertainties. A linear fit 
gives the following mean relation between the new ages 
and those in the GCS: AgeNew = 0.16-1-0.90 x Ageccs (aU 
ages in Gyr). Thus, on average, the largest ages decrease 
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by ^--^10%. This is similar to the differences seen when us- 
ing different stellar models and has negligible impact on 
the interpretation. 




Age old (Gyr) 

Fig. 18. Ages based on the new metallicity and temper- 
ature calibrations and model corrections vs. the original 
ages GCS, using only single stars with ages better than 
25% in both sets. The dotted hne shows the 1:1 relation. 

The only noteworthy deviations occur on the two 
"branches" that can be seen in Fig. [181 This small group 
of stars is located in the "hook" region in the HR dia- 
gram, where an observed point is matched by two differ- 
ent isochrones, one placing the star on the main-sequence 
turnoff, the other on the early subgiant branch. As ex- 
plained e.g. in J0rgensen & Lindegren (|2005p . this leads 
to a two-peaked G-function structure. Small changes in 
the assumed [Fe/H] and T^ff may then change the relative 
height of the peaks, and thus the most probable age of the 
star. For such double-peaked G-functions, one could use a 
weighted mean of the two maximum values rather than a 
simple fit to the highest peak, but the improvement would 
be cosmetic rather than real. 

Finally, we compare our results to two recent sets of 
ages, both computed from the spectroscopic results by 
VF05. One set is given in that paper itself, using a rela- 
tively crude method, while the more recent one by Takeda 
et al. (|2007p was computed with the same Bayesian tech- 
nique as used in the GCS, based on a dense grid of stel- 
lar evolution tracks computed with the Yale-Yonsei code 
(Demarque et al. I2004p . In order to avoid confusion by 
including stars with highly uncertain ages (see, e.g. Fig. 8 
of Reid et al. I2007p . we select only single stars with ages 
better than 25% as given by both sources. 

Fig. [in] shows the remarkable result of this compar- 
ison. While the VF05 ages are on average ~10% lower 




Age VK05 (Gyr) 



Age Takeda (Gyr) 



Fig. 19. Left: New ages vs. ages from VF05. The sample 
consists of single GCS stars with both new and VF05 ages 
better than 25%; the 1:1 relation is shown. Right: The 
same comparison with ages from Takeda et al. (|2007p . 



than ours, those by Takeda et al. (|2007p are in essentially 
perfect agreement with ours. Considering that these re- 
sults are based on totally different observational data and 
analysis techniques, calibrations, and stellar models, the 
agreement is extraordinarily close and inspires a solid con- 
fidence in age determinations from isochrones. Note that 
the differences between Takeda et al. (.2007,) and VF05 are 
due to different computational techniques applied to the 
same data, while the differences between the old and new 
GCS ages are due to improved calibrations, the computa- 
tional method remaining the same. That the systematic 
differences seen in both cases are only of the order of ~10% 
gives further confidence in the robustness of the method. 

For future reference we recall that, temperature and 
metallicity calibrations aside, the age estimates computed 
by HOG were based on an unrealistically small value for 
the uncertainty of Toff. Further, the computations ignored 
the uncertainty in My as well as the temperature offset of 
the models, the a-enhancement of the metal-poor stars, 
and the presence of binaries in the sample. Moreover, 
the bright limit of My=2 excluded many young stars for 
which good ages can be readily determined (see GCS Fig. 
11a), while at the same time including low-mass giants 
for which meaningful ages cannot be derived. Finally, the 
uncertainty of the resulting ages was not discussed. 

8. Results and sample characteristics 

Based on the new calibrations and model corrections dis- 
cussed above, we have computed new Tcft, [Fe/H], dis- 
tances, ages and age uncertainties, and space motions for 
all the stars in the GCS. Table[T]gives the results; the full 
version is available in electronic form at the CDS. We have 
not recomputed Galactic orbital elements for the stars, as 
uncertainties in the assumed (smooth, axisymmetric) po- 
tential are more important than minor revisions of the 
space motions. In the following, we discuss the implica- 
tions of the new data for our understanding of the evolu- 
tion of the Milky Way disk. 

In such discussions, it is crucial to not only employ the 
best possible calibrations from observed to astrophysical 



14 



J. Holmberg et al.: The Geneva-Copenhagen Survey of the Solar neighbourhood II 



Table 1. Sample listing of the recomputed parameters for the first five stars in the GCS catalogue. Ageiow and Age^p 
are the lower and upper 1-a confidence limits on the computed age, respectively. The full table is available in electronic 
form from the CDS (see reference on title page). 
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parameters, but also to understand to what extent the 
criteria used to select the stellar sample may influence 
the conclusions drawn. This is especially important when 
studying subsamples selected on the basis of having one or 
more of the derived parameters available, perhaps within a 
certain level of precision, rather than the entire GCS. Age 
is the most striking example: The samples of stars having 
ages better than 25%, or "well-defined" ages in the GCS 
sense, or any ages at all, are all very different from the full 
GCS sample in very non-random and non-intuitive ways. 

Much effort was spent checking these issues in prepa- 
ration for the GCS. Only overall results were mentioned 
in the paper itself, details being left to the present paper. 
These tests are the subject of the following sections. 

8.1. Simulating the sample selection effects 

To estimate the interplay between the sample selection 
criteria and any errors in the determination of the derived 
parameters in the GCS, we have performed an extensive 
set of numerical simulations based on a synthetic cata- 
logue. The synthetic catalogue assumes a constant den- 
sity of stars in a spherical volume of ~250 pc radius, a 
mixture of 90% thin-disk and 10% thick-disk stars with 
an even distribution in stellar age and evolutionary stage. 
Several distributions in metallicity and kinematics were 
imposed on the thin disk in order to explore the possi- 
ble systematic effects; for the thick disk we assumed an 
asymmetric drift in V of 65 km s— 1 and {au,crv,o'w) = 
(72,43,36) km s-\ 

The star density was normalised such that, after ap- 
plying the appropriate cutoffs, the resulting sample would 
contain about 15 000 stars total, with about 97% thin-disk 
and 3% thick-disk stars, as observed. Each component has 
a distribution in metallicity and kinematics similar to the 
observed values, and the distribution of stellar parame- 
ters and measurement errors is also like the real cata- 
logue. In order to better test the age estimation process 
at faster evolutionary stages, the stars are more uniformly 
distributed along the isochrones than what would result 
from a standard IMF. 

By applying the various calibrations, computations, 
and parameter cuts to the simulated catalogue in parallel 
with the observed sample, we can ascertain in quantitative 
terms which systematic effects may be introduced at each 



stage and whether they have any significant influence on 
the conclusions. 
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Fig. 20. Average volume occupied by single stars with 
"well-determined" ages in the real GCS (full line) and in 
the simulated catalogue with a fixed apparent magnitude 
limit of 8 (dotted). 

As a first example, we show how the volume of space 
sampled by the GCS varies with age. Because older stars 
are, on average, fainter (and cooler) than younger stars, 
and the GCS sample is limited by apparent magnitude, 
the volume sampled by the GCS decreases with stellar age. 
GCS Fig. 24 showed this effect as a function of colour; here 
we illustrate the dependence on age directly by imposing a 
simple apparent-magnitude cutoff at V= 8 and computing 
the average volume occupied by the stars in successive 1- 
Gyr age bins, normalised to that occupied by the youngest 
stars (=100%). 

Fig. [20] compares the result for the simulated sample 
(dotted curve) with the same computation for the real 
GCS (full line). Only single stars with "well-determined" 
ages are included. The GCS adopts fainter limiting mag- 
nitudes for the redder stars in order to achieve volume 
completeness to 40 pc, so the fraction of older stars is 
higher than for the simplified simulation. 

However, the origin of the preponderance of young 
stars in, e.g.. Figs. [23l and [29l is clear: They are included 
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from a volume ^^5 times larger than the stars older than 
~3 Gyr. Accounting properly for these sampling volume 
differences is, of course, crucial when deriving the star for- 
mation history of the solar neighbourhood from the data 
- a study we do not undertake here. 

More specific simulations have been carried out to test 
the most significant conclusions of the GCS; they are de- 
scribed in the following as appropriate in each context. 

9. The metallicity distribution function 

Models of star formation and chemical evolution in the 
galactic disk make predictions of the distribution of heavy- 
element abundances in stars that (could) have survived 
from all stages of the evolution. The classic failure of 
'closed-box' models - the so-called "G-dwarf problem - 
refers to the observed lack of those metal-poor dwarf stars 
that should have accompanied the high-mass stars that 
produced the heavy elements we do observe in the younger 
stars, assuming a constant IMF. One may ask whether the 
metallicity calibration or the way the data are compared 
with the models may cause a spurious difference. 




fMe/Hl T,„ 



Fig. 21. Left: Metallicity distribution function for all stars 
using the new calibration (solid) as well as the original 
GCS values (dotted) . Right: The distribution of Tcfr with 
the new and old calibrations (same symbols). 

Fig. [21] compares the distributions of the GCS stars 
in [Fe/H] and T^s as derived with the old and new cali- 
brations. As seen, the new calibration leads to even fewer 
metal-poor stars in the sample than the old one, but the 
difference is marginal. 

Galactic evolution models typically predict average 
properties of the stellar population in a vertical column 
through the disk at the position of the Sun, while obser- 
vations typically - as in the GCS - refer to a roughly 
spherical volume centred on the Sun. Knowing the mo- 
tions of the stars enables us to allow for this difference. 

We have performed the correction from a volume- 
complete to a column-complete sample in a rigorous man- 
ner, using the self-consistent mass models of the disk de- 
rived by Holmberg & Flynn (j20001 12004} latest version in 
Flynn et al. l2006|) . With this model, we derive the relation 
between volume and column density for each velocity dis- 
persion. Fig. [22] shows the vertical velocity dispersion of 
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[Fe/H] [Fe/H] 

Fig. 22. Left: The W velocity dispersion as a function of 
[Fe/H] for the volume-complete sample (equal-size bins). 
Right: Metallicity distribution function for stars in the 
volume-complete sample (solid); in the column, using the 
fit in the left-hand panel (dashed); and in the closed-box 
model of Casuso & Beckman (|20041 thin line) . 

the GCS stars as a function of metallicity, together with a 
smooth polynomial fit. This fit is then used to transform 
the volume densities to corresponding column densities as 
functions of [Fe/H]. 

Fig. [22| compares the volume and column metallicity 
distribution functions with each other and with a closed- 
box model without instantaneous recycling approximation 
(Casuso & Beckman ,2004.) . convolved with a Gaussian of 
cr= 0.08 to account for the observational error in [Fe/H]. 
Clearly, neither the calibration nor the volume/column 
correction can account for the observed dramatic defi- 
ciency of metal-poor stars relative to the model. 

We cannot comment on the reasons why HOG reached a 
different conclusion from the same data (albeit with a dif- 
ferent metallicity calibration), since the closed-box model 
discussed by H06 was not described there. 

10. The age-metallicity diagram 

The relationship between average age and metallicity for 
stars in the solar neighbourhood ~ the so-called age- 
metallicity relation (AMR) - is probably the most pop- 
ular diagnostic diagram for comparing galactic evolution 
models with the real Milky May. There is, however, no 
consensus on its shape or interpretation. 

The discussion centres essentially on (i) the presence 
or absence of a general slope of the AMR (a gradual in- 
crease in mean [Fe/H] with time), and (ii) how much of 
the scatter in [Fe/H] at any given age is accounted for by 
observational errors and how much reflects the complex- 
ity in the evolution of a real galaxy as compared to most 
current models. It is thus crucial to ascertain whether the 
distribution of stars in the age-metallicity diagram (AMD) 
reflects primarily the intrinsic properties of the sample or 
primarily artefacts of the measurement and parameter cal- 
culation procedures. 

This can be investigated by means of our simulated 
sample. First, we check whether the new calibrations have 
by themselves caused any signiflcant change in the AMD, 



16 



J. Holmberg et al.: The Geneva-Copenhagen Survey of the Solar neighbourhood II 




0.5 - 



10 

Age GCS (Gyr) 




10 

Age New (Gyr) 



Fig. 23. Top: AMR from the original GCS, for single 
stars with ages better than 25% and with E(b-y) < 0.02 
or d < 40 pc. Large dots are mean values in bins with 
equal numbers of stars. Bottom: Same, using the improved 
stellar parameters. 
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Fig. 24. Same as Fig. [2 
sample with d < 40 pc. 
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iU but only for the volume-limited 



using the stars with the very best ages {a{Age) <25%). 
In order to eliminate the bright, distant early F stars in 
the GCS whose high [Fe/H] we suspected to be due to 
overestimated reddening corrections, we have limited the 
sample to single stars within 40 pc or with E(b-y) < 0.02 
mag. As seen, the revised calibrations cause hardly any 
difference in the AMD, neither for the full sample (Fig. 

nor for the strictly volume- limited one (Fig. . 

Perhaps the most obvious structure in the observed 
AMD is the marked increase in mean metallicity and the 
absence of metal-poor stars for ages below 2 Gyr (see Fig. 
[23] and GCS Figs. 27-28). It is important to understand if 
this is due, at least in part, to the selection criteria used 
to define the catalogue sample, or whether it is a genuine 
property of the solar neighbourhood. In the latter case, it 
could be interpreted as a cut-off in the resupply of fresh 
low-metallicity gas, followed by closed-box evolution. 

Fig. mi shows the results of computing the "observed" 
AMR from a simulated sample with a flat input AMR 
(first panel) . The second panel shows the effect of recom- 
puting the ages and metallicities from the synthetic data, 



using the new calibrations and restricting the sample to 
stars with ages better than 25%. This illustrates the vary- 
ing difficulty of determining precise ages over the HR dia- 
gram, especially in the region where isochrones overlap at 
4-5 Gyrs. 

The third panel shows the AMR after applying the 
same blue colour cutoff {b-y > 0.205) and apparent mag- 
nitude limit as used to define the GCS sample. It exhibits 
the same preponderance of young metal-rich stars as the 
observed diagrams (Figs.[23l- [24|) . Thus, the apparent lack 
of young metal-poor stars in the AMD is caused simply by 
the blue colour cutoff of the GCS. None of these selection 
effects seems to have been considered by H06 or by Reid 
et al. pOOT)) . 

In order to verify whether our calibrations or age com- 
putations could introduce (or remove) trends at higher 
ages in the AMD, we imposed tight AMRs (width 0.1 dex 
at all ages) on a subset of the simulated catalogue. In or- 
der to test all combinations of age and metallicity, we im- 
posed both a linear increase as well as a linear decrease of 
[Fe/H] with time. We then computed synthetic colours for 
the stars with random errors as observed, and processed 
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Fig. 26. Alternative input and "observed" AMRs for 
the synthetic catalogue (Fig. [25] showing a flat AMR). 
Left: Input AMR: A 0.10 dex wide band of slope - 
0.036 dex/Gyr. Right: Same, but with the opposite slope 
(+0.036 dex/Gyr). 
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Fig. 25. AMR simulations for the synthetic sample. Top 
left: Input AMR for a volume complete synthetic sam- 
ple, flat within each of the thin and thick disks (dots and 
open circles, respectively). Top right: "Observed" AMD 
for the synthetic sample, using the new calibrations with 
the synthetic observations and retaining only stars with 
resulting ages better than 25%. Bottom left: "Observed" 
AMD as above, but after imposing a blue colour cutoff (6- 
y > 0.205) and apparent magnitude hmit as in the GGS. 
Bottom right: Derived ages vs. the "true" input values. 



Subtracting the estimated observational error in each case 
(0.10 and 0.07 dex) yields an intrinsic ("cosmic") scatter 
in [Fc/H] at a given age of 0.20 dex - identical to the 
result by Edvardsson et al. (|1993|) from high-resolution 
spectroscopy. 

We conclude the discussion by reiterating the impor- 
tance of the sample selection on the resulting AMR. If cer- 
tain types of stars are excluded a priori, the gross shape 
of the AMR may be predetermined whatever new metal- 
licities or ages may be observed. Fig. [27] illustrates this, 
using two popular sources, VF05 and Edvardsson et al. 
p-993 ). Both were selected for specific purposes, and both 
are prone to selection effects affecting the AMR. 

Edvardsson et al. selected evolved F-type dwarfs in or- 
der to be able to determine isochrone ages, excluding a 
priori any old, metal-rich and therefore redder (i.e. G) 
stars that might have populated the upper right-hand part 
of the AMR. VF05 studied stars used in planet searches 
and thusavoided stars with weak lines, i.e. hot metal-poor 
stars in the opposite corner of the AMR. These facts are 
clearly discussed in the papers, but make these and other 
inherently biased samples unsuitable for discussions of the 
general AMR of the solar neighbourhood. 



these synthetic stars with the new calibrations in the same 
manner as in the GGS, including the colour cutoffs, and 
retaining only stars with ages better than 25%. 

Fig. [26] compares the input and recovered AMD for 
both of these cases. As can be seen, the slope of the in- 
put AMR is faithfully reproduced in both cases, despite 
the inevitable scatter introduced by the observational er- 
rors. We thus conclude that the absence of a significant 
mean slope of the data points in Fig. [24] reflects the true 
situation of the solar neighbourhood. 

Focusing on the stars older than 4 Gyr in Fig. [24l 
we find a mean [Fe/H]= -0.24 dex and standard devia- 
tion cr= 0.22 dex for the original GGS data (top panel, 
70 stars), and a mean [Fe/H]= -0.21 dex and a= 0.21 
dex using the new calibrations (bottom panel, 75 stars). 



10.1. AMR slope vs. a radial metallicity gradient 

Recently, Rocha-Pinto et al. (|2006p interpreted an ob- 
served variation of mean metallicity with the difference 
between the mean orbital radius of the star (Rm) and that 
of the Sun (R©) as evidence for a tight age- metallicity 
relation in the thin disk. However, because their values 
of R,n range from 6 to 9 kpc, many of their stars must 
have large velocities and orbital eccentricities, usually as- 
sociated with thick-disk stars. This was already noted by 
Edvardsson et al. (|1993p . although they did not use the 
explicit term "thick disk" . 

A much more natural explanation of the apparent ra- 
dial variation in metallicity is that the observed stellar 
sample is a mix of thin-disk stars dominating close to the 
solar circle, and thick-disk stars dominating at both high 
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Fig. 27. AMR for two pre-selected samples using different 
metallicity and age determination methods. All ages have 
errors below 25%. Top left: The VF05 sample, using our 
data. Top right: The Edvardsson et al. sample, using our 
data. Bottom left: The VF05 sample, using their data. 
Bottom right: Same, using Takeda et al. ages. 
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Fig. 28. Radial variation in [Fe/H] for the single GCS 
stars (dots; new calibration) , and for a simulated catalogue 
with the AMD of Fig. [25] (line), i.e. with no correlation 
between age and metallicity in either the thin or the thick 
disk, but a radial metallicity gradient of -0.09 dex/kpc 
in the thin disk. We obtain a similar overall relation as 
Rocha-Pinto et al. I|2006p . but without any AMR in either 
thin or thick disk. 



and especially low Rm- To verify this, we have compared 
the observed radial metallicity distribution for single stars 
in the GCS (using the new calibrations) with that of our 
simulated catalogue, which has the flat age-metallicity dis- 
tribution of Fig. [25l i.e. no correlation between age and 
metallicity in either the thin or the thick disk. However, we 
have added a radial metallicity gradient of -0.09dex/kpc 
for the thin-disk stars, close to the observed values (see, 
e.g., discussion in Sect. 6.2 of the GCS). Fig. [28] shows 
that this reproduces the observed radial metallicity vari- 
ation without any variation of metallicity with age. 

We have deliberately made no attempt to fit the de- 
tailed structure in the Rm - [Fe/H] distribution. This de- 
pends on the complex structure of the distribution of the 
U and V velocities (see Figs. [29] - [30] and GCS Fig. 20), 
which exhibits a multitude of kinematic groups not ac- 
counted for in a simple diffusion picture. As an example, 
the peak in [Fe/H] at Rm — 6.8 kpc is associated with 
the C Herculis stream discussed in the GCS and again by 
Famaey et al. ([2005]) and Bensby et al. ([2007]) . 
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11. Age-velocity relation and disk heating 

As demonstrated earlier in this paper, our new calibra- 
tions cause no significant systematic changes in the ve- 
locities and ages of the GCS stars (apart from the overall 
~10% age scale reduction). Thus, no major revision of the 
GCS results on the evolution of the kinematics of the disk 
with age is to be expected. However, the size of our sample 



Fig. 29. U, V and W velocity vs. age for the 4065 single 
GCS stars with a{Age) < 25%). 
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Fig. 30. U-V diagrams for the GCS subsamplc of Fig. 
separated into four age groups. 



Fig. 32. W velocity distributions for the GCS subsample 
of Fig. [29l separated into four age groups. No significant 
substructure is seen. 
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Fig. 31. V-W diagrams for the GCS subsample of Fig 
separated into four age groups. 



allows us to examine that evolution in greater detail than 
was done in the GCS. 

For this task, we select the cleanest sample of GCS 
stars, i.e. the single stars with a{Age) < 25% and with 
complete space velocity data. Note that, due to our im- 
proved age determination, there are now 4065 stars in this 
class as compared to 2852 in the original GCS ~ an increase 
by almost 50%. 



Fig. [29] shows the observed space velocity components 
as functions of age, while Figs. [30] and [31] show them in 
the U-V and V-W planes, separated into four groups by 
age. Like Figs. 20 and 30 of the GCS, they illustrate the 
significant substructure in the U and V velocity distri- 
butions that persists over a wide range of ages (see also 
Famaey et al. 120051 and the study of the Hercules stream 
by Bensby et al. l2007|) . In contrast, the W velocities show 
no deviation from a random distribution in any age group 
(Fig. [32]), suggesting that different heating mechanisms 
are at work in the plane of the disk and perpendicular to 
it. However, it must also be noted that phase mixing is 
stronger for vertical motions due to the Galactic poten- 
tial, and the orbital periods shorter, resulting in a more 
efficient smoothing of structure. 

A detailed discussion of the mechanisms affecting disk 
star orbits in and perpendicular to the plane is beyond 
the scope of the present paper, but we can consider two 
basic questions, i.e. (i) whether the velocity dispersion of 
disk stars continues to increase during the lifetime of the 
disk or whether it shows signs of a plateau or saturation in 
some age interval, and (ii) what functional form appears 
appropriate for the rising parts of the age- velocity relation 
(AVR). For this discussion, we divide the sample by age 
into 30 bins with equal numbers of stars (i.e. 135 in each 
bin), against the only 10 bins used in the GCS; this allows 
us to follow the evolution in greater detail. 

Fig. [33] shows the resulting AVR, which shows a 
smooth, general increase of the velocity dispersion with 
time in both U, V, and W. Fitting power laws while ex- 
cluding the three youngest and three oldest bins, we find 
exponents of 0.38, 0.38, 0.54 and 0.40 for the U, V, W and 
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Fig. 33. Velocity dispersions vs. age for the GCS subsample of Fig. [29] (4065 single stars with a{Age) < 25%; new 
calibrations). The 30 bins have equal numbers of stars; the lines show fitted power laws. The 3 youngest and oldest 
bins have been excluded from the fit. 



total velocity dispersions - slightly larger than the values 
derived in the GCS, as expected from the change in age 
scale. 

This, however, is clearly not the full story of the phe- 
nomenon known as "disk heating" . If that term is taken to 
refer to an increase in the random motions of disk stars, 
it is clear from Fig. [30] that the present distribution of 
the U,V velocities is not the result of pure "heating"; 
non-random processes are at work as well. Several possi- 
ble mechanisms have been described in modern literature, 
based on detailed simulations, but all result in velocity 
dispersion increases with time that can be approximated 
by a power law. This is true, e.g. for the simulations that 
use transient spirals as the heating agent (De Simone, Wu 
& Tremaine l2004p or produce velocity structure and heat- 
ing of the disk by interaction between two spiral systems 
(Minchev & Quillen i2006p or between a spiral and the bar 
(Chakrabarty]2007l). 

To identify any underlying "pure" heating mechanism 
within the disk, one would have to remove the major non- 
random substructures, such as the Hercules stream or the 
SiriusUMa, Coma, and Hyades- Pleiades branches (GCS; 
Famaey et al. 120051 Bensby et al. I2007p . This, however, 
introduces a certain arbitrariness as to the choice of stars 
to be removed; and it appears more interesting to us to 
identify the cause(s) of the substructures that dominate 
the U,V plane than to quantify any remaining minor ef- 
fects. 

On the other hand, the observed heating of the verti- 
cal (W) velocities, which show no signature of non-random 
mechanisms, is more efficient than found in simulations. 



E.g. Hanninen & Flynn (|2002p find an exponent of only 
0.26 using molecular clouds and need massive black holes 
to reproduce the observed heating rate, in agreement with 
several earlier studies. The black-hole hypothesis is incon- 
sistent with other observational data, but the extra ver- 
tical heating might be due e.g. to infalling satellites with 
dark matter substructure (Benson et al. 2004p. In any case, 
further theoretical work on the heating mechanisms affect- 
ing the vertical velocities in the disk will have to satisfy 
the observational constraint illustrated in Fig. [33] 

This conclusion contrasts with that by Quillen & 
Garnett (|200ip . who claimed that the heating of thin-disk 
stars saturates at a maximum velocity dispersion ~3 Gyr 
after the birth of a star; the oldest stars were assumed to 
belong to the thick disk and heated by a different mecha- 
nism. This led Freeman & Bland-Hawthorn (|2002p to sug- 
gest that remnants of early mergers might still be traced 
as dynamical groups in the disk today (see Helmi et al. 
120061 for an actual application). 

Fig. [33] shows no such saturation in our data, but the 
levelling-off of the velocity dispersion in the two oldest bins 
might be taken as weak evidence for a constant velocity 
dispersion in the thick disk. When judging the merits of 
the two results, it should be remembered that the sample 
discussed here is some 20 times larger than that discussed 
by Quillen & Garnett (|200ip and the sampling errors cor- 
respondingly smaller. 

Finally, we want to verify whether our age computation 
techniques might bias the slopes seen in the age-velocity 
relation (AVR). To do so, we have imposed a fixed AVR on 
our synthetic catalogue, assuming purely Gaussian veloc- 
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ity distributions and adopting the slopes of cr([/, V, W) vs. 
time as derived above. The thick disk contribution in the 
simulation is at the same level as in the observed sample, 
about 3% after applying the apparent-magnitude cutoff. 
We then recomputed the ages and age errors for all stars 
in the synthetic sample, using our new calibrations and re- 
taining only stars with (new) ages better than 25%, sorted 
the stars into 30 equal-size bins, and recomputed the ve- 
locity dispersion in each bin. 

The input and reconstructed AVRs are compared in 
Fig. [34l Symbols with error bars show the "observed" 
points with recomputed ages and age bins. As can be seen, 
all the reconstructed velocity dispersions are quite consis- 
tent with the input AVRs over the range 1.5-8 Gyr. There 
is thus no evidence for any error or bias being imprinted 
on the AVR from the age estimation process. 
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Fig. 34. Age- velocity relations for the synthetic sample. 
Solid lines: Adopted input relations for thin-disk stars. 
Symbols and error bars: Recomputed ages and velocity 
dispersions in 30 equal-size bins for stars with ages better 
than 25%. Dotted lines: Fitted power laws. 



12. Thick disk vs. thin disk 

In the literature, different criteria are used to distinguish 
thick-disk stars from the more numerous thin-disk stars 
in the Solar neighbourhood, based on kinematics, detailed 
chemical composition, or age. There is, however, no real 
consensus on what limiting values to adopt in any of these 
parameters, and the allocation of any individual star to ei- 
ther disk component is often ambiguous. Questions of cur- 
rent interest include the fraction of thick-disk stars in the 
Solar neighbourhood as well as the existence or otherwise 
of an AMR in the thick disk. 



The thick disk of the synthetic sample is well visible 
in the reconstructed AMR and AVRs, especially among 
the oldest stars, where it dominates. In order to estimate 
the fraction of thick-disk stars in the GCS, we calculated 
the fraction of thick-disk stars in the simulated catalogue 
with V-velocities below -100 km/s and [Fe/H] > -1. This 
sample contained 27% of the total number of thick-disk 
stars in the catalogue, and no thin-disk stars were found 
in this velocity-metallicity domain. In the real, observed 
GCS we find 67 stars in this region. Thus, if we assume the 
same thick-disk velocity distribution as in the simulation, 
we find a total of 248 thick-stars; as the total number of 
single stars with complete space velocities in the GCS is 
8479, the thick-disk fraction is thus 2.9%. 

As regards a possible AMR in the thick disk we note 
that, interestingly, the "observed" AMR for the simulated 
thick-disk stars in Fig. [^5] shows a pronounced AMR of 
slope -0.022 dex/Gyr, although the input "true" AMR had 
a constant metallicity between the assumed thick-disk age 
limits of 11 and 12 Gyr. This effect would have been even 
more pronounced if our strict quality criteria on accepted 
ages had not been applied (as in H06), or if the thick- 
disk sample had contained a small fraction of stars with 
younger ages and/or higher metallicity. 

Samples of field thick-disk stars can be constructed 
by combining different membership criteria. For nearby 
field stars, a kinematic criterion is most commonly used, 
based on a decomposition of the local velocity distribution 
into Gaussians corresponding to the thin and thick disk 
and assuming values for the asymmetric drift and U, V, W 
velocity dispersions for each. Membership probabilities for 
either disk are then assigned for each star, based on the 
observed space motion. 

First, as shown in Fig [30] and GCS Fig. 20, the actual 
local velocity distribution is nothing like a two-Gaussian 
model; in particular, the C Herculis stream occupies a po- 
sition near the overlap between the thin and thick disks. 
Second, Fig. [34] shows that the practice of assuming a 
single set of velocity dispersions for the thin disk, regard- 
less of age, ignores the effects of the continuing dynamical 
heating of the disk. As general thin-disk samples are typ- 
ically dominated by relatively young stars (see Figs. [23l 
and [29]) . the mean values of a{U,V,W) will tend to be 
underestimated, enhancing the probability that thin-disk 
outliers will be classified as thick-disk stars. 

If the thick disk is indeed very old, as is most com- 
monly assumed, the kinematic parameters characterising 
the thin disk should be chosen to correspond to the old- 
est thin-disk stars (see Fig. [34]). It is interesting to note 
that Bensby et al. (|2004p . who identify thick-disk stars 
with intermediate ages and chemical properties, assume 
significantly lower values of (y{U, V) for the thin disk than 
Reddy et al. (|2006p . who find a smaller fraction of such 
stars and consider them to be outliers from the thin disk. 

Given the inherent ambiguity of kinematic population 
indicators, the enhanced [a/Fe] ratio observed in thick- 
disk stars may be the best membership criterion for indi- 
vidual field stars. In order to study the AMD of the thick 
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disk from the cleanest possible sample of stars, we there- 
fore selected all stars in common between the GCS and 
those kinematically classified thick-disk stars from Reddy 
et al. 11)06) that also have [Ti/Fe] > 0.25. 

Fig. [551 shows the AMD for this "clean" thick-disk sam- 
ple, using our new ages. We see a tight grouping of stars 
with ages around 11 Gyr and [Fe/H] ^ -0.55, and no trace 
of a slope. Taken together. Figs . [25l and [35l put into ques- 
tion the claims of a significant age-metallicity relation in 
the thick disk by, e.g. H06 or Bensby et al. PDM]) . 

Note that a strong, spurious correlation of age with 
metallicity would result if the Padova isochrones had not 
been corrected for the increasing discrepancy between the 
original isochrones and the observed unevolved stars at 
decreasing metalhcity (see Sect. [5]). This is also the case 
if the a-enhancement of metal-poor stars is ignored when 
choosing the Z parameter of the models. 



0.5 - 



-0.5 
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Fig. 35. AMD for thick-disk stars with [Ti/Fe] > 0.25 
from Reddy et al. (|2006p . Only single stars with ages bet- 
ter than 25% and with E(b-y) < 0.02 or d < 40 pc are 
shown. 



13. Conclusions 

We have redetermined the basic calibrations used to infer 
astrophysical parameters for the GCS stars from uvby pho- 
tometry. The result is a substantial improvement in the 
determination of effective temperatures, especially for the 
hotter stars, and of absolute magnitudes (i.e. distances), 
as well as minor corrections to [Fe/H]. The GCS redden- 
ing calibration appears to be essentially correct. These 
new results are obtained by drawing on the large body 
of 2MASS K photometry, high-resolution spectroscopic 
abundance analyses, and accurate Hipparcos parallaxes 
that are now available. 

Using the improved astrophysical parameters, we have 
recomputed the ages and age error estimates for the GCS 
sample (Tabled]). In the process, we have compared re- 
sults from a variety of stellar models, finding no large dif- 



ferences between models for typical GCS stars. We have 
also recomputed the temperature corrections needed for 
the models to agree with the observed unevolved main 
sequence at metallicities below solar. 

The resulting ages correlate well with those published 
in the GCS, but are on average ~10% lower. The recent 
ages by VF05 for ~ 400 stars in common are some 10% 
lower still, while those by Takeda et al. (j2007p are essen- 
tially in perfect agreement with ours, although they were 
also based on the spectroscopic temperatures and metal- 
licities by VF05. This independent verification places the 
determination of isochrone ages on a very firm basis, pro- 
vided adequate precautions are taken. 

We note that, because of a small offset in the GCS uvby 
photometry for the Hyades stars and the non-standard He 
abundance of this cluster, the Hyades cannot be used to 
check metallicity or age scales for field stars such as those 
in the GCS, as assumed by H06. The revised results given 
here for Hyades and Coma stars are based on standard 
photometry whenever available. 

The revised [Fe/H] values change the observed metal- 
licity distribution only marginally; it is still in strong dis- 
agreement with the prediction of closed-box galactic evo- 
lution models. 

In preparation for the discussion of Galactic relations 
involving stellar ages, notably the age-metallicity and age- 
velocity diagrams, we have performed extensive simula- 
tions of the effects of our selection and computation pro- 
cedures by applying them to synthetic catalogues with 
properties closely resembling those of the GCS, but with 
a variety of specified intrinsic properties, such as the AMR 
and AVR. We find that our methods faithfully recover the 
input relations within the observational errors, without 
introducing spurious trends or other systematic effects. 

The observed AMR retains the general features of that 
of the GCS, i.e. little or no variation in mean metallicity 
with age in the thin disk, plus an admixture of perhaps 3% 
thick-disk stars, and with a large and real scatter in [Fe/H] 
at all ages. We find no evidence for a significant AMR in 
the thick disk. Note that these conclusions are only ro- 
bust when backed by careful end-to-end simulations of the 
properties of the sample and associated selection criteria 
(notably the limits by apparent magnitude and the blue 
colour cutoff), and by meticulous attention to the details 
of the computation of the stellar ages and their uncertain- 
ties. Investigations using simpler approaches (e.g. H06 or 
Reid et al. I2007P may well reach different conclusions. 

A significant result of the GCS was the demonstra- 
tion that the dynamical heating of the thin disk continues 
throughout its life. We confirm this from our revised data 
set, and with substantially higher time resolution than in 
the original GCS. Our simulations also confirm that no 
bias has been introduced in the slope of the AVR for the 
individual U, V, W velocity components by our age com- 
putation technique. 

We conclude that kinematic parameters taken to rep- 
resent the thin disk in thick-disk/thin-disk separations 
based on observed velocities should be chosen to corre- 
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spond to the oldest thin-disk stars. Using average values 
for thin-disk samples dominated by stars much younger 
than the thick disk may lead to contamination of the thick- 
disk samples by kinematic outliers from the thin disk, es- 
pecially when such non-Gaussian kinematic features as the 
Hercules stream are present. 
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